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Abstract 

We solved the Frenkel-Kontorova model with the potential V(u) = — ^\X\(u — 

Int[u] — |) 2 exactly. For given |A|, there exists a positive integer q c such that 

for almost all values of the tensile force a, the winding number uj of the 

ground state configuration is a rational number in the g c -th level Farey tree. 

For fixed cu = p/q, there is a critical A c when a first order phase transition 

occurs. This phase transition can be understood as the dissociation of a large 

molecule into two smaller ones in a manner dictated by the Farey tree. A 

kind of "commensurate-incommensurate" transition occurs at critical values 

of a when two sizes of molecules co-exist. "Soliton" in the usual sense does 

not exist but induces a transformation of one size of molecules into the other. 
PACS number(s):03.20.+i, 05.45.+b, 64.60.-i 
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I. INTRODUCTION 



Frenkel-Kontorova (FK) model describe a system of particles moving in an infinite se- 
quence of potential wells. In the limit of shallow wells, the particles are kept at an equilibrium 
distance by a tensile force a. Such models are also widely studied in the context of circle 
maps when certain periodicity exists in the locations of the potential wells |T] . A stationary 
configuration in the FK model corresponds to an orbit in the circle map. In the study of 
FK models, however, one is particularly concerned with the minimum energy configurations. 
The ground state is a "recurrent" minimum energy configuration that is characterized 
by a winding number u, the inverse of which, 1/u, gives the average number of particles per 
well. Hence a ground state with rational u = p/q Q consists of molecules that are composed 
of q particles and p wells. 

It is well-known that two types of phase transitions occur quite generally in FK models 
||. The first type is the commensurate- incommensurate transition which occurs at critical 
values of a when the enthalpy for the creation of solitons or anti-solitons vanishes 0. The 
second type occurs at critical values of A characterizing the strength of the potential wells 
and corresponds, in the language of circle map, to the breaking of an invariant circle into a 
cantorus. In the ground state with an irrational winding number u, the allowed positions of 
the particles in the potential wells change from being the entire period in the "unpinned" 
phase to being a cantor set in the "pinned" phase. Correspondingly, the hull function f u (x), 
from which one obtains the position u n of the n-th particle by u n = f u (nuj), turns from a 
smooth function into a no-where-differentiable function P,|^,|g-|8| . 

These two types of phase transitions, however, are not the only phase transitions that 
could occur in a FK model. In particular, when the potential possess an "internal struc- 
ture", e.g., having multiple local minima in a period, one will encounter first order phase 
transitions at critical values of the height of a local minimum [§,|nj. At the critical point, 
the filling fraction of the corresponding well (the local minimum) changes from a value sub- 
commensurate to the winding number u to a value non-sub-commensurate to uj |IIJ]. As in 
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the case of commensurate-incommensurate transition, this new type of first order transition 
is nucleated by solitons ||. 

In this letter, we would like to present yet a new type of first order phase transition by 
solving exactly an FK model with a potential which is concave almost everywhere. The 
phase transition is not nucleated by solitons. Rather, it can be understood as the breaking 
up of a large molecule into two smaller ones, in a manner dictated by the Farey tree. 

The stationary configuration of an FK model is determined by the equation 

u n+1 - 2u n + u n -i = XV' (u n ), (1) 

where u n denotes the coordinate of the n-th particle. The model we consider has the potential 

1. / , r , l^ 2 



K(u) = ~A(u-Int[ii]--J , -4<A<0, (2) 
where Int[w] equals the largest integer not larger than u. The case of A > had been 



thoroughly investigated before by several authors [[11] . 

The recurrent stationary configurations are characterized by a winding number uj and 
can be specified by a hull function /^(x) so that 

u n = fu{nu + a), (3) 

where a is a constant phase. The hull function satisfies + 1) = f u (x) + 1. 

For u — p/q, the hull function is a step function with q steps in a period. Depending on 

whether f u (x) is continuous at integer values of x, there are two solutions: 

9-1 

fw{x) = u n(q) (lnt[x + nu] - lnt[nu]j , (4) 



n=0 



and 



/«(*) = [/„(x + ^)+/ w (x-l)]/2,, (5) 

where u n (q) = tan | esc ^ cos(n — |)%, x — arccos(l — |A|/2) and lnt[x] = (Int[x + ] + 
Int[x _ ])/2. The configuration described by f w (x) has one particle sitting at the potential 
bottom in a given period. fu,(x) can also be expressed as 
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fu( x ) = J2 u n(q) [lnt[x + nu + — ] - Int[no; + — ]J , (6) 

where v n = [y n + u n i)/2 with n' = n + g x (mod g), < n' < g — 1, and gi is determined 
by the Farey tree in the following way. Given the irreducible fraction u = p/q, ui\ = P\/q\ 



and co>2 = is the unique irreducible pair such that []12 

p = pi + P2, q = qi + 52, ^1 < u; < ^2 (7) 

and 

pgi = pig + 1, p 2 g = 2*72 + 1, P2gi = Pig2 + l- (8) 

We shall call the "g-th level Farey tree" to be the Farey tree where irreducible fractions 
with denominator larger than q are truncated. The average energy per particle corresponding 
to the stationary configurations described by f u (x) and fu(x) are respectively 

= y~t E "M {I - (™ ~ \ ~ intM) 2 } (9) 

and 

f (u) = - ^^(g). (10) 

8g 

One notices that at x — Xg = ^/^j ^0(9) vanishes and the two configurations become 
degenerate. Eq. (|T0|) also indicates that for < x < Xg; the ground state configuration with 
winding number u = p/q is given by fu,(x). 

It is instructive to learn how the configuration described by f w (x) evolves as x increases 
from zero. Let u n = /^(nw — l/2g) for < n < g — 1. As x increases from zero to Xqi the 
Un's change from v = v\ = • • • = = 1/g to a first zero value occurring at z/ . When 
this happens, the two particles Uq and u qi touch the bottom of the potential wells. As x 
increases further, these two particles are pinned at the bottom and we can no longer use 
fu){%), which ceases to be an increasing function, to describe this configuration. Instead, the 
gi particles from uq to u qi -i are described by an u f Ul section", i.e., a consecutive gi particles 



chain with the first particle sitting at the potential minimum in the f UJl configuration, while 
the q 2 particles from u qi to u q _\ are described by an f W2 section. For a general to — p/q, if we 
regard an f u section as a "molecule" of size (q,p), i.e., composed of q particles and p wells, 
then Xq can be seen to be a critical point of a first order transition when a molecule of size 
(q,p) is just about to break up into two molecules of sizes (qi,pi) and {q 2 ,p 2 ) respectively. 
There three type of molecules, whose corresponding u>i, u, and u 2 are related as consecutive 
Farey fractions, co-exist at the critical point. 

For x > Xq-> fu( x ) ceases to describe the ground state configuration with winding number 
uj. One can show from Eqs. (§) and (|10D that 

qtVfa) + g 2 f (u; 2 ) - q^(u) = ± tan g cot ^ cot ^ cot ^ (11) 

o Z Z Z Z 

and the right hand side is negative when x > Xq- The molecules of size (q,p) dissociate and 
the ground state now is a mixture of / Wl and f W2 sections with the right proportion. 

When x continues to increase and reaches x q = it/ max(gi, q 2 ) = min(x gi , x 92 ), say 
Xq — Xqn the molecule of size (gi,Pi) starts to dissociate into even smaller molecules of sizes 
(q 2 ,P2) and (q x - q 2 ,Pi -p 2 )- Note that (pi —p 2 )/ (?i - q 2 ), Pi/qi, and p 2 /q 2 are consecutive 
Farey fractions in the (g — l)-th level Farey tree. The process continues until all the particles 
are located at the bottom of the potential well. This occurs at x > 7r /2- 

It is interesting to observe that for x q < X < Xq- even though f^x) no longer describes 
the ground state configuration, it remains an increasing function until x reaches x q . Hence 
molecules of size (q,p) may be regarded as an unstable resonance state of two smaller 
molecules of sizes (<?i,pi) and (q 2 ,p 2 ) when x hes in the interval (XqiXq)- 

In summary, for given x, there exists a positive integer q c such that it/ (g c +l) < X < ^/q c - 
The ground state configurations can only be composed of the f w sections with u in the q c - 
th level Farey tree in the following way. For an arbitrary irrational winding number uj 
or rational but not in the g c -th level Farey tree, we can find a unique pair of consecutive 
fractions U\ and uo 2 in the g c -th level Farey tree such that uj\ < uo < uo 2 . The ground state 
configuration with this given u can be constructed with a fraction fi of particles associated 
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with the f Ul sections and a fraction f 2 associated with / Wa sections such that 

fl + h = 1, fl^l + 1*2^2 = ^- (12) 

Moreover, the average energy per particle of this ground state configuration is given by 

*M = ^^l) + ^^ 2 ). (13) 
Co>2 — CUi dJ 2 — ^1 

Adding or removing one particle from the system will turn p\ / W2 sections into p 2 
sections or vice versa. This can be seen from Eq. (|8|) which implies that pxP2(l/^i — 
I/1V2) = 1 • Starting from pure f Ul sections, i.e., a f Wl configuration, we can approach 
a fu> 2 configuration by adding particles one by one. Solitons and anti-solitons in the usual 
sense of local "defects" do not exists. Effort to create them merely induces transitions 
between f wi and f W2 sections. If we insist on calling "defect" an f Wl section in a background 



of f W2 sections, this "defect" will carry a "fractional charge' 
It follows from Eq. (|TT| ) that 



P2 ' 



) y( Ul ) + [ -5- ] f > f(w) for X < X q , (14) 

which indicates that ^ e (u;) is a convex function of cu. Indeed, for given Xq c +i < X < Xg c > this 
can be shown by making use of Eq. ( |13|) and proving the convexity of ^f(u) for any three 
consecutive fractions uii < uj < u 2 in the g c -th level Farey tree. It follows from Eq. @ that 
Pi + Vi = Kp and gi + g 2 = K( ?o f° r some positive integer k. Similar to Eq. ( pi]) one can 
derive that 

?i*(o;i) + <? 2 *(^2) - Kq^{uj ) = - tan ^ cot cot cot (15) 

o z z z z 

The right hand side is positive and this completes our proof. 

Up to now, we have been discussing the ground state configuration for given uj as \ is 
varied. Introducing the tensile force term —<joj into Eq. ([13]), we obtain the enthalpy of 
the system. By minimizing this enthalpy with respect to u, we obtain the phase diagram 
in the x ~ a plane, as shown in Fig. 1. One can see that there exist tricritical points at 
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Xq = 7r/q, for each q, where three types of molecules co-exist. For x < Xgj the lowest 
enthalpy configuration is locked to uo = p/q for 

gO^-j^i) <(j< V(uo 2 ) - 

UO — UO\ ~ ~ L0 2 — 00 

where cu 1 < uo < uo 2 are consecutive fractions in the g c -th level Farey tree with Xq c +i < X ^ 
Xq c - When x — Xqi ^ ne s ^ e P shrinks to a point and we have 

(7 = = , (17) 

UO — U>i L0 2 — LO 

which gives the equations for the tricritical point. In general, for two consecutive Farey 
fractions ooi and oo 2 , cr = ^^5^^ gives the equation for the co-existent curve of the two 
phases corresponding to molecules of sizes (qi,Pi) and (q 2 ,P2)- 

In conclusion, we have shown that the solvable FK model considered in this work illus- 
trates some important physics accompanying a first order phase transition. In particular, we 
demonstrate explicitly how the Farey tree dictates the structure of these phase transitions. 
Even though the model is one-dimensional and the potential is not smooth, the underlying 
physics could very well be realized in nature. 
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